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We investigate the scissors modes in binary mixtures of degenerate dilute quantum gases, for 
both Bose-Bose and Bose-Fermi mixtures. For the latter we consider both the superfluid and normal 
hydrodynamic and collisionless regimes. We analyze the dependence of the frequencies of the scissors 
modes and their character as a function of the Bose-Fermi coupling and the trap geometry. We show 
that the scissors mode can reveal a clear trace of the hydrodynamic behavior of the Fermi gas. 



I. INTRODUCTION 

In recent years, the improvement of the trapping tech- 
niques has allowed for the creation of multi-component 
Bose-Einstein condensates (BECs), formed by atoms in 
different internal (electronic) states 0, The multi- 
component BEC, far from being a trivial extension of 
the single-component one, presents novel and fundamen- 
tally different scenarios for its ground-state wavefunction 
an d excitations [1,013 In particular, it has been 
experimentally observed that the BEC can reach an equi- 
librium state characterized bythe phase separation of the 
species in different domains p|- 

During the last few years, the mixtures of fermions 
and bosons, from now on called Fermi-Bose (FB) mix- 
tures, have also attracted a growing attention. Although 
the initial interest has been mostly motivated by the 
possibility to achieve sympathetic cooling towards the 
Bardeen-Cooper-Schriffer (BCS) transition H, [^}, re- 
cently, the rich physics of FB mixtures has become it- 
self one of the central topics of the physics of ultracold 
gases. Various properties of the FB mixtures have been 
analyzed, including the phase separation between bosons 
and fermions the existence of novel types of 

collective modes 0, Tl3t llil Il5| , the appearance of ef- 
fective Fermi-Fermi interactions mediated by the bosons 
El H E 113, the collapse of the Fermi cloud in the 
presence of attractive interactions between bosons and 
fermions [3, the physics of ID FB mixtures or FB 
mixtures in optical lattices [2(i|. 

The analysis of collective excitations is an excellent 
tool to analyze the effects of interactions in ultracold di- 
lute gases. Zero-temperature excitation frequencies have 
been extensively studied in the case of Bose-Einstein con- 
densates (BECs) in dilute alkali gases up to a high degree 
of experimental accuracy |2ll |22| . These experimental 
results are in excellent agreement with the theoretical 
predictions obtained from first principles p3l l24l |25| . 

Among the different collective modes the so-called scis- 
sors mode has attracted considerable attention. This 
mode is well known in nuclear physics |26j . where it cor- 
responds to the out-of-phase rotation of the neutron and 



proton clouds. In the context of ultracold atomic gases, 
the scissors mode is achieved after a sudden tilt of the 
trap, and it consists in the harmonic motion of the atomic 
cloud around the new trap axis. This particular collective 
excitation, predicted in Ref. 1271 . has been successfully 
observed experimentally j2^, |29| |. 

The present paper is devoted to the analysis of the 
scissors modes in mixtures of degenerate atomic gases, 
both for the case of a multicomponent BEC, and for a 
FB mixture. In the latter case, we analyze both when the 
Fermi gas is in hydrodynamic and in collisionless regime. 
We discuss in particular the variation of the frequencies of 
the scissors modes (as well as their associated geometry) 
as a function of the Bose-Bose (or Bose-Fermi) coupling 
constant. 

The paper is organized as follows. In Sec.|n]we analyze 
the case of a Bose-Bose mixture. Sec. 11111 is devoted to 
the analysis of the scissors modes in FB mixtures, first for 
the case of a hydrodynamic Fermi gas, and later when the 
Fermi gas is in the normal phase. We finalize in Sec. IIVI 
with our conclusions, and the discussion of the results. In 
App.^we have included some analytical results, which 
have been employed in different calculations throughout 
the paper. 



II. BINARY BOSE-EINSTEIN CONDENSATE 

In this section we consider the case of a BEC formed by 
two different components {1,2}, confined in a harmonic 
external potential 



Vi(r) 



iX 2 + w 2 , 



2 i 2 

y +u zv , 



(i) 



where i — 1,2. We are interested in the scissors modes of 
this particular system after a sudden tilt of the trap axis. 
We assume that the chemical potential is larger than the 
energy separation between the trap levels (Thomas-Fermi 
regime). Therefore the effects of the quantum pressure 
are negligible, and the system can be analyzed using the 
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hydrodynamic equations for superfluids 



^+V-(n lVl ) = 



(2) 

(3) 
i^- + v(^y-^ + yi+.gii?n+ 5 i 2 n 2 ) =0(4) 

= 0(5) 



— + V • (n 2 v 2 ) 
9vi 



9v 2 /m 2 2 

m 2~Qf + V + ^2 + 01211 + 322n 2 



Since we are interested in the oscillations of the 
cloud in the x-y plane after a tilting from the equilib- 
rium solution, it is convenient to adopt the coordinates 
x\ = xcos(0j(i)) — y sm(0i(i))) and y[ = xsin(9i(t)) — 
y cos(8i(t))). We consider small angles 9i, and con- 
sequently we can perform the Ansatz m(t,x,y, z) — 



nig(xi,yi,z) ~ n g (x,y,z) + 6i(yd x 



xd v 



c(x,y,z) 



where rij is the density, the velocity field, to, the mass, 
and Vi the external potential, associated with the i-th 
component. The coefficient gij = 2ith 2 ! ay / 'rhij is the 
coupling constant associated with the interatomic inter- 



Note that the time dependence is entirely contained in 
the parameters 9\ and 9 2 (see Fig. 1). 

After multiplying the continuity equations by xy and 
integrating, we obtain the equation for the asymmetry of 
the cloud of the i-th component 



dt {xyh = 



yv x )i 



(8) 



actions between the components i and j, where my and Note that the asymmetry para meter is related to the an- 



are the corresponding reduced mass and scattering gular d i S p lacem ent 9 t in the form {xy), = (y 



length, respectively 




Figure 1: Schematic description of the two angles of rotation 
involved in the scissors modes. 



Neglecting the quantum pressure term in a single com- 
ponent BEC demands that Na/l^o ^> 1, where N is 
the total number of particles and Iho = y/H/miD, where 
u) = Lo x Lo y ijj z . This approximation breaks down at the 
border of the Thomas-Fermi density profile, where the 
density is small. In the case of a two-component BEC, 
one must additionally require that the two species are 
not strongly phase separated. Otherwise, at the domain 
wall, the kinetic energy dominates the mean-field poten- 
tial, and therefore the quantum pressure is certainly not 
negligible. 

Imposing a stationary solution without flow (v = 0), 
we obtain 

Vi(r) + gnn lg (v) + g 12 n 2g (r) = fn (6) 
Va(r) + gwnig(r) + g 22 n 2g (r) = ^ 2 , (7) 

from where we can easily calculate the ground-state den- 
sity profiles n ig (r). 



where (. . . ) g denotes the expected value in the equilib- 
rium state. Neglecting terms of second order in the an- 
gular displacements, we obtain 



yvx) = J (y,x,o) 



dv 

fit 



nd 3 r. (9) 



Substituting into Eq. Q the corresponding Euler equa- 
tion for the i-th components, and combining with Eqs. 
©, we obtain, in the linear regime 



BuOi — -An 0i 
-S 22 6* 2 = A 2 i9\ 



Ai 2 2 
A 22 2 , 



(10) 
(11) 



where the coefficients Ay and Bi can be found in Add.IAI 
By imposing a solution of the form exp{— iuoi), we find 
the frequencies of the different scissors modes 



AuB 2 + A 22 5i ± y/(A n B 2 + A 22 5i) 2 - dot [A] 

B\B 2 

(12) 

where we have introduced the matrix A = {Ay}. The 
corresponding eigenvectors {9f , 9^} provide information 
about the nature of the modes. The two components 
oscillate in-phase when 9\9 2 > 0, and in counter-phase, 
when 9\9 2 < 0. 

In Fig. 2 we depict our results for a particular set of 
parameters (see figure caption). Both in-phase and out- 
of-phase modes depend on the particular value of g\ 2 . 
Notice, however, that the in-phase mode is less affected 
by the change of g\ 2l since the interaction energy be- 
tween both gases is just slightly modified during the mo- 
tion. In particular, if the trapping frequencies for both 
components are the same, the in-phase mode remains in- 
dependent of the interaction coupling constant g\ 2 . The 
counter-phase mode is in any case strongly affected. 

The in-phase mode is always the lowest energy state 
for g\ 2 < 0, since it maximizes the overlapping. For 
.912 > 0, the out-of-phase mode becomes the state with 
lowest energy, as long as the two components present a 
maximum in the trap center, since in that case the out-of- 
phase mode minimizes the overlapping. For a given value 
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.912 = 3i2' the density profile of the component with low- 
est central density (say the component 2) acquires zero 
curvature at the trap center, therefore, for g 12 > 9%, one 
of the components presents a minimum at the trap cen- 
ter, and the in-phase mode becomes once more the lowest 
energy state. For w x , y , Z ;i = u x ,y,z;2> ffS = 9n, whereas 
from the Thomas-Fermi solutions © and 10, one can 
easily obtain that in general, if (3 = oj x>y>z .x/ui XtVtZ -2, 
512/511 = 2 {m 1 /m2)- This dependence of the modes 
has been confirmed by our numerical simulations (see 
Fig. EJ). 

III. BOSE-FERMI MIXTURE 

In this section we analyze the scissors modes for FB 
mixtures. Two different regimes are considered: (i) The 
case in which the fermions are in the hydrodynamic 
regime, and (ii) the case in which the fermions are in the 
normal phase in the collisionless regime. Concerning the 
case (i) , we would like to note that the scissors modes for 
a FB mixture in the normal phase in the hydrodynamic 
regime, coincide with those modes of the superfluid hy- 
drodynamics. In Sec. IV we comment on this issue in 
more detail. 



A. Hydrodynamic regime 
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Figure 2: (a) Frequencies of the scissors modes for a binary 
BEC as a function of the intercomponent coupling constant 
g\2 for Ni = N2 = 10 6 , ui x ,y,z;2 = 0-97u x>y , zl i, equal masses 
m, uiy-i = 1.5wa;;i, and uj z -i = 2oi x; i. The dimensionless inter- 
action energies are gu/HuJx;ilx-i = 1 and <?22 = 2gn, where 
i'ji — h/7TUJ x -i. The frequencies are indicated in units of 
— (u)x-,i + UJ 'ya) 1 ^ 2 ■ (b) Character of the modes, provided 
by £ = 2#i#2 /{o\ + 62), as a function of g\2. The in-phase 
(out-of-phase) mode is characterized by £ > (£ < 0). 



In a two-component Fermi gas, if the typical collisional 
time becomes much larger than the trapping period, the 
Fermi gas is said to enter into the hydrodynamic regime. 
As for the Bose-Bose mixture, the FB mixture in the 
hydrodynamic regime can be well described by means of 
the hydrodynamic equations 



dt 

dvy 
dt 




— + V • (n 6 v b ) = 

V • (n f v f ) = 
Vb + gtrib + gbffif 

nib 

Vf + g b fn b h 2 



(13) 
(14) 
(15) 



TO/ 



2m 2 f 



(67r 2 n/ 



,2/3 _ 



v/A(VAv/). 



(16) 



Note that in our calculations only one set of equations 
is considered for the Fermi gas. This is justified if the 
difference between the trapping frequencies and the con- 
centrations of both Fermi components is negligible |30| . 
Note also that we have employed the chemical potential 
for a free Fermi gas, ~ n /^ 3 ) since in the quantum degen- 
eracy regime due to Pauli-blocking (away from Feshbach 
resonances) the corrections due to Fermi-Fermi interac- 
tions are small. 

As for the Bose-Bose mixture, we obtain the ground 
state density profile in the Thomas- Fermi approximation, 



co/co, 



Fl 
1.025 



\ 




Figure 3: (a) Frequencies of the scissors modes for a FB mix- 
ture, where the Fermi gas is in the hydrodynamic regime, 
in units of uip± = {(jJ^.j + ujy.j) 1 ^ 2 , as a function of the 
intercomponent coupling constant gbf, (b) Character of the 
modes, provided by £ = 26b0f/(9i + (??), as a function of gbf- 
The in-phase (out-of-phase) mode is characterized by £ > 
(£ < 0). In the figures we consider the case Nb = Nf — 10 6 , 
equal masses m, U) x>y>zi b = 1.05lj IiHiZ; /, lu x .j — 0.1u z -j, 
ujy-f = 0.995lj z; /, and gb/hu) z jl 3 z .j = 1. 
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and introduce two time dependent parameters 9 b and 6 f , 
which describe the small angular displacements of the 
bosonic and fermionic cloud, respectively. From the hy- 
drodynamic equations we obtain the equations of motion 
for the angles 



C 1 6 b = D n 9 b + D 12 9 f 
C 2 f = D 21 9 b + D 22 6 f . 



(17) 
(18) 



Although the Eas. l|17[) and (|18|) are formally the same as 
Eqs. (|10fl and l|ll|) . the coefficients are obviously differ- 
ent. The form of these coefficients is detailed in ApplAl 
The frequencies of the scissors modes are obtained in a 
similar way as in Eqs. i|12|l . A typical result is depicted 
in Fig. [3J Note that similarly to the case of Bose-Bose 
mixtures, and due to similar reasons, the out-of-phase 
mode becomes the less energetic state for < g b f < g^j- 
Note that, contrary to the Bose-Bose case, Eq. <|16f) 
allows for the possibility of rotational flow, which could 
be present if the Fermi gas is in the hydrodynamic regime, 
but still in the normal phase. However, since we are only 
interested in small deviations with respect to equilibrium, 
the rotational term, as well as the kinetic energy term, 
are neglected, since they provide corrections of second 
order in the angular displacements. 



d)/<0_ 





Sbf^b 



Figure 4: (a) Frequencies of the scissors modes for a FB mix- 
ture, where the Fermi gas is in the collisionless regime, , in 
units of ljf± = (wj.j + ojy.j) 1 ' 2 , as a function of the inter- 
component coupling constant / ; (b) Character of the modes, 
provided by £ = 29 b 9 f /(9l + 6)), as a function of gbf . The in- 
phase (out-of-phase) mode is characterized by £ > (£ < 0). 
Same parameters as in Figure [3] are considered. 



in the collisionless regime. In this case the bosons still fol- 
low the superfluid hydrodynamic equations (|13f) and 11511 . 
but the fermionicgas follows the so-called Boltzmann- 
Vlasov equation 



dt 



1 



-V P E/> • V v / = 0, 



(19) 



where / is the phase-space distribution function for the 
fermions. The external potential, Uf(r) = Vf(r) + 
g b fn b (r), contains both the harmonic confinement for the 
fermions (V/) and the mean field term due to the inter- 
action with the bosons. The fermionic density is defined 
as rif(t,r) = J d 3 vf(t,r,v). 

Combining the hydrodynamic equations for the bosons 
and the Boltzmann-Vlasov equation for the fermions, we 
obtain the following set of equations 



dt 
d_ 

di 



(xy) b = (xv y + yv x ) b 
(xv v +yv x ) b = 



(xd y + yd x )(^v! + V b 



g b n b + g b fn f 



— (a;?/)/ = {xv y ) f + (yv x ) f 

^(xvy)f = (v x v y ) f + -L(x[d y (V f + g bf n b )]) f 

-jl(v v x}f = (vxVy)f + —{y[d x {Vf+g bf n b )]} f 
at mf 

^( v xV y )f = ([(v x dy + v y d x )(V f +g bf n b )]) f , 



(20) 

(21) 
(22) 
(23) 

(24) 

(25) 



where V b is the harmonic confinement for the 
bosons, (x) b = $ x{r)n b {t, r)d 3 r/N b , and (\) f = 
J x(r,v)/(i,r,v)d 3 rd 3 v/iV / . 

After a rather tedious algebraic manipulation we arrive 
to the following equations of motion for the displacement 
angles (see App. 1X1 : 

L x e b = M 11 6 b + M 12 6f (26) 

L 2 9 f = M 21 9 b + M 22 6 f + M 23 r] (27) 

r? = M u e b + M a2 9 f +M 3 3V- (28) 

Note that in the collisionless regime a third variable 
couples with the displacement angles, namely r\ = 
J v x v y fd' i rd 3 v/Nf, which is the average value of the de- 
formation in momentum space [2?|. Note that in the 
hydrodynamic regime rj = 0, since in the hydrodynamic 
case the momentum distribution remains always spher- 
ically symmetric. A typical case is depicted in Fig. 21 
Note the appearance of three scissors modes, and the 
mixing of the bosonic and the collisional fermionic modes. 



B. Collisionless regime 



IV. CONCLUSIONS 



If the typical time for Fermi-Fermi collisions becomes 
much smaller that the trapping period, the system enters 



In this paper we have analyzed the scissors modes ob- 
tained after a sudden tilting of the trap axis for the case 
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of binary mixtures of quantum gases. We have studied 
the scissors modes for the case of a BEC with two compo- 
nents, and especially the dependence of the modes with 
the intercomponcnt coupling constant. Whereas for at- 
tractive intercomponcnt interactions the in-phase mode 
always presents the lowest energy, we have shown that 
for the case of repulsive interactions the character of the 
modes depends on the density profile of the ground state, 
and in particular on the overlapping of both components. 
If both components present the maximal density at the 
trap center, the lowest energy corresponds to the out-of- 
phase mode. At a given critical coupling constant one of 
the components develops a minima at the center, and the 
in-phase mode becomes the less energetic. A similar rea- 
son explains the behavior of other collective excitations 
in binary mixtures |15| . 

In the second part of the paper we have focused our 
attention on the case of a mixture of degenerated bosonic 
and fermionic gases. We have first analyzed the case in 
which the fermionic gas is in the hydrodynamic regime, 
in which the scissors modes present a similar behavior as 
for the case of binary condensates. We have completed 
our analysis with the detailed study of the case in which 
the Fermi gas is in the normal collisionless phase. In that 
case, the system presents three different scissors modes, 
due to the asymmetry of the velocity distribution of the 
fermionic component. 

Finally, we would like to stress that the difference be- 
tween the superfluid hydrodynamic regime and the colli- 
sional hydrodynamic regime is provided, from the macro- 
scopic point of view, by the irrotational constraint for the 
velocity field in the superfluid case (V A v = 0), related 
to the quantization of the flux and the appearance of 
quantized vortices. On the contrary, this condition is not 
necessarily fulfilled in the collisional hydrodynamics. As 
discussed above, the scissors mode is excited after a sud- 
den change of the trap axis. In the absence of friction 
such a sudden change cannot induce a rotational flow. 
Consequently, the scissors modes in collisional hydrody- 
namics do not differ from those expected in the superfluid 
case. The differences between collisional and superfluid 
hydrodynamics in FB mixtures can be revealed extend- 
ing to the FB mixtures similar ideas as those discussed 
for Fermi gases in Ref. jl^l. This problem will be the 
subject of further studies. 

Note added After the completion of this work, we have 
learned of a related work performed by Kasamatsu et al. 
|33|| on the problem of scissors modes in binary conden- 
sates, in which the authors employ different techniques as 
those discussed in this paper. In Ref. [33| special empha- 
sis is paid to the nonlinear coupling between quadrupole 
and scissors modes. 



Acknowledgments 



gases, ESF PESC BEC2000+, and the Ministero 
dell'Istruzione, deU'Universita e della Ricerca (MIUR). 
M. R. acknowledges the hospitality of the Hannover 
group. L. S. and P. P. wish to thank the Alexander 
von Humboldt Foundation, the Federal Ministry of 
Education and Research and the ZIP Programme of 
the German Government. M. R. and P. T. thank the 
Academy of Finland (Project Nos. 53903, 48445) for 
support. 

Appendix A 

The coefficients of Eqs. (f 101) and can be ob- 

tained after a somewhat tedious calculation following the 
method described in the text, and acquire the form: 



1 
1 



Bi = — I xyn' lg d 3 r 



B 2 = jj~ I xyn 2g d r 



(Al) 
(A2) 



A u = J n lg (yd x +xd y )n' lg d 3 r (A3) 



Au = H j n lg (yd x +xd y )n' 2g d 3 r (A4) 

A 2 i = || J n 2g (yd x +xd y )n' lg d 3 r (A5) 

A22 = ^Jn2 9 (yd x +xd y )n' 2g d 3 r, (A6) 

where rn g denotes the equilibrium density, and n' ig — 
(yd x - xd y )n ig . 

Similarly, we can obtain the coefficients of Eqs. i|17|) 
and 02 



1 



Ci = — / xyn' lg d 3 r 



C 2 = jr Jxyn 2g d 3 r 



N f 
9b 



(A7) 
(A8) 



D n = j?- J n bg (yd x + xd y )n bg d 3 r (A9) 
D 12 = 9 jj-J n bg (yd x + xd y )n' fg d 3 r (A10) 



9bf 



D21 = ^ / n fg (yd x + xd y )n' bg d 3 r (All) 



iff 
N f 



1-hi = ^- / n fg (yd x +xd y )-—[^d\ (A12) 



2ri< 



fg 



where g// = fr 2 (67r) 2 / 3 /2m/. 

Finally, let us briefly discuss the derivation of 
Eqs. (|2"rJj) . (|2"T|) and J5SJ). We assume a velocity field for 
the fermions of the form 



We acknowledge support from Deutsche Forschungs- 
gemeinschaft (SFB 407), the RTN Cold Quantum 



v x (t,r) = aiy 
v y (t, r) = a 2 x 



(A13) 
(A14) 
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where the time dependence is contained in the parame- 
ters a 1 and a 2 . We additionally define 

(v x v y ) = T] (A15) 
(xv y + yv x ) = (3. (A16) 

It is possible to rewrite Eqs. I|2() [l -l|25 jl in the form 



Ei9b 


= a 


(A17) 


$ 


= F\9 b + F 2 9f 


(A18) 


(G 2 -G 1 )9 f 


= G±ai + G 2 a 2 


(A19) 


Gidi 


= rj + Q^b + Q^f 


(A20) 


G 2 a 2 


= v + Q2iO b + Q 22 9 f 


(A21) 


V 


— + R 2 a 2 


(A22) 



where 



£1 
Fi 
F 2 
Gi 
G 2 

Q12 



1 

N~b 
9b_ 



xyn' bg d 3 r 

n-bgiydx + xd v )n bg d 3 r 



jjr J n bg (yd x + xd y )n' fg d 
y 2 n fg d 3 r 
x 2 nf g d 3 r 



1 
W f 
1 



u x 1 . 9bf 1 „ 



(A23) 
(A24) 
(A25) 
(A26) 
(A27) 
d 3 r (A28) 



Qn = jf f J n Sg yd x n' bg d 3 i 



Q 



22 



W V I , 9bf 1 

N- f Xynf9 + Nf~ nf t 1 ynbg 



= -rr- I n fn xd. 



Ri = u 



2 



L fg XU y U bg t 



(A29) 



d 3 r (A30) 



(A31) 



R 2 = 



Gx + jf f n bg yd v n bg d 3 r (A32) 

(A33) 



lo1G 2 + 



7 
9bf 
N f 



n b gxd x nbgd 3 r 



The coefficients of Eqs. (|26I28|) are of the form 



Li = 


Ei 




L 2 = 


G 2 - 


Gi 


Ma = 


Fi 




M 22 = 


F 2 




M 2 i = 


Qn- 


hQ 2 i 


M 22 = 


Q12- 


h<?22 


M 23 = 


2 




M31 = 


RiQ 
Gi 


11 + 


M 32 = 


RiQ 
Gi 


12 + 


M 33 = 


Ei 

Gi ^ 


R 2 
' G 2 



R 2 Q 2 i 

G 2 
R 2 Q 22 
G 2 



(A34) 
(A35) 
(A36) 
(A37) 
(A38) 
(A39) 
(A40) 

(A41) 
(A42) 
(A43) 
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